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The naturally weak spin-orbit coupling in Graphene can be largely enhanced by adatom deposi¬ 
tion (e.g. Weeks et al. Phys. Rev. X 1, 021001 (2011)). However, the dynamics of the adatoms also 
induces a coupling between phonons and the electron spin. Using group theory and a tight-binding 
model, we systematically investigate the coupling between the low-energy in-plane phonons and the 
electron spin in single-layer graphene uniformly decorated with heavy adatoms. Our results provide 
the foundation for future investigations of spin transport and superconductivity in this system. In 
order to quantify the effect of the coupling to the lattice on the electronic spin dynamics, we com¬ 
pute the spin-flip rate of electrons and holes. We show that the latter exhibits a strong dependence 
on the quasi-particle energy and system temperature. 


I. INTRODUCTION 

Two-dimensional topological insulators exhibiting the 
quantum spin Hall effect (QSHE) hold great potential 
for applications in low-power consumption electronics.' 
However, so far very few such materials are available. 4,5 
This makes the quest for new materials exhibiting the 
QSHE a very important research enterprise. Early in the 
history of the field, it was pointed out by Kane and Mele 1 
that graphene, (l ’ if endowed with a sizable spin-orbit 
coupling (SOC), would exhibit a robust QSHE. However, 
the strength of the SOC in graphene is just a few tens of 
micro-volts, 8 ’ 9 ’ 12 ’ 1 which makes the QSHE in graphene 
extremely fragile and, in practice, inaccessible in actual 
experiments. 

Nevertheless, after its isolation in 2004, graphene, a 
monolayer of carbon atoms arranged in a honeycomb lat¬ 
tice, has become an extremely popular 2D material be¬ 
cause it can be cheaply produced in large-size samples 
displaying large carrier mobilities. Furthermore, the car¬ 
rier charge and density can be easily tuned by gating. In 
recently years, several methods for enhancing the SOC in 
graphene have been put forward. ’ 1-16 They largely rely 
on adatom deposition, and can lead to either predomi¬ 
nantly Rashba-type SOC, 1 predominantly intrinsic-type 
SOC, 1,1 or both, ■ depending on the adatom type. Cas¬ 
tro Neto and Guinea pointed out that chemisorption of 
atomic species like Hydrogen or Fluorine, which bind to 
the carbon atoms by inducing sp 2 -sp 3 hybridization, lo¬ 
cally induces a spin orbit coupling of the Rashba-type, 
which can be as large as several tens of meV. This pro¬ 
posal has been recently verified experimentally by Bal- 
akrishnan and coworkers, 1 who measured the spin-Hall 
effect (SHE) in lightly hydrogenated graphene devices. 

In addition, it was proposed by Weeks et al. that, 
by deposition of certain heavy adatoms like Indium and 
Thallium, the strength of the intrinsic-type SOC can be 
largely enhanced in single-layer graphene. The enhance¬ 
ment is due to second-neighbor hopping that is medi¬ 
ated by the spin-orbit-splitted p-orbitals of the adatoms, 


which strongly hybridize with the unoccupied 7r-levels of 
graphene. Such an enhancement of the intrinsic SOC can 
lead to a band gap of the order of several tens of meV, and 
therefore, to a much more robust QSHE. Furthermore, it 
was noted that, for these atomic species, the adatoms 
find their lowest absorption energies at the center of the 
hexagons on the honeycomb lattice (the so-called H po¬ 
sition) . Indeed, first principle calculations by Weeks and 
coworkers 1 found that the heavy adatoms at the H po¬ 
sition have binding energies of —0.525 eV, being the en¬ 
ergy difference between H and B position (i.e. the bridge 
between two neighbor carbon atoms) 80 meV. In this 
configuration, the Rashba-type SOC, which is detrimen¬ 
tal for QSHE 1 is entirely absent at the K and K' points 
(but grows linearly with the crystal-momentum separa¬ 
tion to those points). Deposition of other atomic species, 
such as transition metal elements like Iridium and Os¬ 
mium, has been also studied. In this case, the details 
of microscopic mechanism by which SOC is induced on 
the graphene layer by transition-metal atoms containing 
d orbitals is different from the atomic species studied by 
Weeks and coworkers, 11 ’ for which the p-orbitals play the 
dominant role. 

It is worth pointing out that recent experiments 1 1 car¬ 
ried both in CVD graphene (where Copper atoms exist in 
the form of a residue resulting from the fabrication pro¬ 
cess) and in intentionally decorated graphene with noble 
metal clusters of e.g. Cu, and Au, support the idea that 
decoration can enhance the SOC in graphene. As in the 
case of Hydrogen, 1 decoration produces a robust SHE 
in standard Hall-bar devices made with CVD and exfo¬ 
liated graphene samples that have been decorated with 
adatoms. Indeed, single-layer graphene is particularly 
prone to the existence of resonances in the neighborhood 
of the Dirac point induced by adatoms (see 13 and refer¬ 
ences therein). This has been recently shown to lead to 
a sizable enhancement the SHE. 19 ’ 20 

Nevertheless, the possibility of a proximity-induced 
SOC in graphene still remains controversial, as a recent 
experiment 21 found no experimental evidence of such 
effect induced in the magnetoresistance, quantum Hall 
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effect and non-local spin Hall effect of graphene dec¬ 
orated with Indium. However, other studies 22 ’ 2 ' 3 have 
found evidence of such proximity induced SOC by Gold 
and Lead in the quasi-particle spectrum measured ei¬ 
ther through angle-resolved phono-emission or scan¬ 
ning tunneling spectroscopy . 2 1 Further studies are there¬ 
fore needed in order to better clarify the complexities of 
the coupling between the adatoms and the electrons on 
the graphene layer. 

One important concern with the decoration of 
graphene by heavy adatoms is that the latter are often 
physisorpbed rather than chemisorpbed. Thus, the bind¬ 
ing of the adatoms to the layer is rather weak and a native 
estimation of the phonon frequencies of the modes associ¬ 
ated with the motion of the adatom yields a relatively low 
characteristic energy of ~ 10 meV . 38 Since the adatom 
is responsible for the SOC, when its position fluctuates, 
it will affect the electron spin. Added to the relatively 
low-energy of the relevant phonons, it can be expected 
that the lattice dynamics in adatom decorated graphene 
can have a substantial effect on the spin-dynamics of the 
low-energy electrons and holes. The precise way in which 
those low-energy phonons couple to the electron spin de¬ 
gree of freedom is the subject of this work. In addition, 
in order to assess the effect of this coupling on the spin- 
dynamics, we have computed the phonon-contribution to 
the spin-flip rate of electrons and holes. 

The rest paper is organized as follows. In Sec. II, we 
use symmetry arguments as well as a simple model to 
obtain the phonon spectrum of decorated graphene. In 
Sec. Ill, we review the most important facts about the 
tight-biding model for the electronic structure of single¬ 
layer graphene decorated with a periodic array of heavy 
atoms sitting on the H position introduced in Ref. [16]. 
In Sec. IV, the coupling of the low-energy phonons to the 
electron spin is obtained both from symmetry arguments 
and the tight-biding model. The consequences of this 
coupling for the electron spin-flip rate are discussed in 
Sec. V. Finally, we discuss our results in Sec. VI and 
summarize our main conclusions in Sec. VII. Some of 
the more technical details of the calculations have been 
relegated to the Appendices. 


II. PHONON SPECTRUM OF ADATOM 
DECORATED GRAPHENE 

In order to obtain the spectral degeneracies of the 
phonon spectrum around special points of the Brillouin 
zone (BZ), we first rely upon group theory. Thus we 
analyze the symmetry properties of graphene covered 
by a periodic layer of adatoms of the type considered 
by Weeks et al." The symmetry group constrains the 
types of phonon-mediated spin-orbit coupling. Since the 
adatoms sit periodically at the H position on top of the 
graphene layer (see Fig. 1, for example), each unit cell 
contains two carbon atoms and one adatom. Note that, 
in this configuration, the mirror symmetry where z —> —z 
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TABLE I: Character table for the vibrational representation 
of the group C" 1V 

is absent. Therefore, the point group of graphene with 
adatoms on the H position is Cq v , which has 12 elements, 
whose equivalence classes are { E , 2Ce, 2 G 3 , C 2 ,3cr a , 3&' a } 
where C m describes a rotation by 27r/m by an axis (z) 
perpendicular to the graphene layer, <r 0 describes reflec¬ 
tions containing the £ axis that leave the A and B sub¬ 
lattices invariant, and cr' a describes reflections that swap 
the two sublattices. 2! 

Generally speaking, two kinds of phonon modes affect 
the dynamics of electrons at low energies. Phonons near 
the r point can scatter low-energy electrons within the 
same valley. In addition, phonons at K(K' ) points can 
also scatter electrons between valleys because K' is equiv¬ 
alent to —K and 3K is a reciprocal lattice vector umklapp 
scattering (i.e. (K) e + (K) ph +G = (K')e' for G = -3K). 
In order to treat the inter- and intra-valley phonon scat¬ 
tering at the same time, it is convenient to follow Ref. [24] 
and triple the unit cell (see Fig. 1), which means that two 
translations t a 1 and t a2 are factored out from translation 
group and must be included in the point group Cq v ■ This 
yields the group Cg v = C 6v + t al C 6v + t a 2 C ev . This ap¬ 
proach maps the K and K' points to the T point of a lat¬ 
tice with a unit cell that is three times larger. Thus, the 
states at K , K' become degenerate and we can deal with 
them simultaneously. Therefore, the possible inter- and 
intra-valley couplings between the phonons and the elec¬ 
tron spin, which in addition must respect time-reversal 
symmetry, can be classified according to the irreducible 
representations ( irreps ) of the Cq v . 

A. Symmetry analysis of the phonon modes 

First, let us study the classification of in-plane phonon 
modes according to the symmetry operations of Cq v . The 
tripled unit cell contains 9 atoms, i.e. 6 carbons and 3 
adatoms. When one atom is mapped onto an equivalent 
atom by a symmetry operation in Cg , it contributes 1 to 
the character of the transformation. For example, under 
G 3 , only the 3 adatoms are invariant, and therefore the 
character for this element equals 3. The character table 
for the different symmetry transformations in the tripled- 
unit cell lattice is shown in Table I. By applying the or¬ 
thogonality theorem, this representation can be reduced 
to 2Ai + B\ + E[ + G'. Taking into account that a two 
dimensional vector like the atom displacement u = (x, y) 
transforms according to the E\ irrep of Cq v , the vec¬ 
tor of in-plane vibrations belongs to an 18-dimensional 
reducible representation. The latter can be decomposed 
into {2A 1 +B 1 +E[+G')xE 1 = 2E 1 +E 2 +E[ +£'+2G', 
where we have used that- Ai x Ei = Ei,B\ x Ei = 
E 2 ,E[ x E 1 = G', G' x Ex =£'+£?'+ G'; G' is a 
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FIG. 1: Enlarged unit cell contains 9 atoms, 6 of which are 
carbon (in red and blue, for A and B sublattices, respectively) 
plus 3 adatoms (in yellow), ai and a 2 are two vectors of unit 
cell. 


four-dimensional irrep and the E \, , E [, E' 2 are two di¬ 

mensional irreps. 

In order to better identify which of the above modes 
belong to the T point of the BZ, we can go back to the 
original three-atom unit cell rather than using the tripled 
unit cell. The original unit cell transforms according to 
the Cq v group. This group is also the little group for 
the T point (the little group for K and K' is C^ v ). The 
representation of the phonon modes at T can be decom¬ 
posed into 2Ei + Ei , which corresponds to two acoustic 
modes (transforming according to one of the Ei irreps) 
and two optical modes transforming according to £d-|--Ev¬ 
idence, it follows that the K and K' points contribute a 
total of 12 optical modes, which transform according to 
E[ + E' 2 + 2 G' in the (reducible) C' ( ' vibrational repre¬ 
sentation. 


B. Central-force model for the in-plane phonons 

We next confirm and extend the previous findings 
based on the symmetry analysis by explicitly obtaining 
the dispersion and wavefunctions for the phonons in the 
entire BZ. To this end, we use a simple model that as¬ 
sumes only nearest-neighbor harmonic forces: ,26 

Hph = 2 ^ A + — 'y ' ■ [ui — Uj)) , (1) 

i ij 


the unit vector along the direction connecting the i and 
j sites, Ui = ( Xi,yi ) is the displacement of an atom at 
site *, and kij = k is spring constant that determines the 
strength of the bond between two neighboring carbons. 
However, kij = k' determines the force between neighbor¬ 
ing adatoms and carbons. The above central-force model 
describes, for small atomic displacements, the changes in 
kinetic and elastic energy. 

Let us next introduce a six-component vector of 
displacements, u(R), whose components are u(R) = 

[xa(-R) I Ia(R) x b (R) y B {R) x 0 (R) y 0 (R)] T , where 
xatVAiXb and ys correspond to the in-plane displace¬ 
ment of the carbon atoms from their equilibrium posi¬ 
tions on the A and B sublattices, respectively, and Xq, yo 
are displacements of adatom from its equilibrium H po¬ 
sition; R = nidi + 7120,2 {ni,ri 2 being integers) spans 
the Bravais lattice (see Fig. 1), whereas Si describes 
the position of the atoms within the unit cell (i.e. the 
basis). Note that we neglect out of plane dynamics, as¬ 
suming that the system is on a substrate, which damps 
out the flexural phonons. Using periodic boundary con¬ 
ditions, the normal frequency in the three-atom unit 
cell can be computed by expanding Ui,\(R) in running 
waves (A = x, y): 


,»(*) = yf £ £ (2) 


fl=l q 


where e^xiq) describes to the pattern of atomic dis¬ 
placements for the pth normal mode and satisfies ej)(<jf) ■ 
£p,'(q) = a six-component vector whose 

components are a(<?))■ The explicit form of Eq. (1) 
in momentum space is given in the Appendix A. In terms 
of second quantization u M (q) = x/%/2ujp{b ll (q)+b^A—q)), 
the phonon Hamiltonian becomes 


H ph = Ew M (q)6j,(q)6 M (q). (3) 


The phonon spectrum along the high symmetry points is 
shown in Fig. 2. Two special symmetry points, T and 
K (or K'), are marked by yellow lines and phonon modes 
at both points labelled according to their irreducible rep¬ 
resentations. 

Let us next discuss the properties of the lowest energy 
phonon modes. The acoustic (Goldstone) modes whose 
energy vanishes at T enter the intra-valley electron- 
phonon coupling only through the space and time deriva¬ 
tive of the mode displacement. 21 This is because an uni¬ 
form translation of the crystal cannot affect the motion 
of electrons (neglecting the coupling to the substrate). 
On the other hand, the lowest-energy optical phonon at 
r transforming according to the E\ irrep can give rise 
to an electron-phonon coupling that is linear in the dis¬ 
placement field u(R). The square of the frequency for 
this mode is given by the expression: 


where atomic mass nii = m for the carbon atoms, 
TUi = m a for the adatoms (we take m a > to), pij is 
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FIG. 2: Phonon spectrum along high symmetry points. When 
the unit cell is tripled, the K and K' points are mapped to the 
r point of the new BZ. The corresponding modes, from lowest 
to highest energy, transfrom according to the irreps G', E' 2 , G' 
and E[, respectively. In addition, the F point of the original 
(three-atom) unit cell contributes with six modes (of the them 
acoustical) transforming according to the Ei,Ei,E 2 irreps. 


Note that, for a heavy adatom (i.e. m a m ), Wp —> 

3k'/2m, which only depends on the coupling between the 
adatom and the carbon atoms, k !, and the carbon atom 
mass, m. Typically, heavy ad-atoms are physisorpbed 
and therefore k' <C k, meaning that the E\ is expected 
to have a rather low-energy, as we have indeed assumed in 
the calculation shown in Fig. 2. The polarization vectors 
for the Ei mode read: 
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For crystal momentum q near K(K'), phonons can 
scatter electrons from one Dirac point to the other. The 
low-lying degenerate optical phonons transforming ac¬ 
cording to the G' irrep have squared frequencies given 
by: 
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The corresponding polarization vectors are: 

/ * \ 


ei (K') = B 


e 2 (K’) = B 


i 
1 

—i 
1 
0 

_ 4 a _ 

\ — b-\- y/b 2 —4ac ) 

-1 

i 

-1 

—i 

4 a 

-b-\-\/b 2 — 4ac 

0 


( 8 ) 


( 9 ) 


where B is the normalization constant, b = —2k'm + 
(k + k')m a ,a = k'^/m a m, and c = —k'^Jm a m. The 
polarization vectors at K are obtained from time-reversal 
symmetry which requires that ei(K') = e{(K), C 2 (A'') = 
Note that, for k k 1 and m a m, ei(iv') —>• 
[0,0,0,0,0,1] T and e 2 {K') [0, 0,0, 0,1, 0] T . That is, 

in the limit of weak adatom-carbon coupling and large 
adatom-carbon mass ratio, the mode is dominated by 
the motion of the adatoms whereas the displacement of 
the carbon atoms becomes negligibly small (cf. Fig 3b) 
When considering the modes at K , K' points it is con¬ 
venient to take their real and imaginary part combina¬ 
tions. The pattern of E\ mode of F point and G' mode of 
K and K' are shown in Fig. 3. For Indium on Graphene, 
the value of adatom-carbon elastic constant has been es¬ 
timated to be 1 k! ss 0.02 Ry/bohr 2 ~ leV/A 2 between 

carbon and Indium atoms with k « 50 eV/A" between 
neighbor carbon atoms. Hence, we obtain Huir — 25 meV 
and hwK — 11 meV by Eq. (4) and (7), which is an energy 
scale comparable to the calculated band gap M = 42 
meV for this system. 


III. ELECTRONIC STRUCTURE OF 
DECORATED GRAPHENE 

A minimal tight-binding model for a single-layer 
graphene sheet uniformly covered by adatoms was given 
in Ref. [16]. In what follows, we briefly review the 
main results obtained in Ref. [16] that are relevant to 
our derivation of the phonon-SO coupling from this mi¬ 
croscopic model. This will also allows us to introduce 
the continuum-limit Hamiltonian that describes the elec¬ 
tronic states at low energies. 

The Hamiltonian of the model contains three terms: 

Ho = H g + H a + H c . (10) 

H g describes the hopping the between the n orbitals of 
the two nearest neighbor carbon atoms: 

H 9 = 4cr'+h-C.). 

(rr') 


\m a m 


( 11 ) 
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FIG. 3: (a) Phonon modes 2E1 + E 2 at F and (b) Phonon 
modes G' = E\ + E 2 at K and K'. In (b) the small displace¬ 
ments of carbon atoms are not shown. 

The term allowing electrons to tunnel between the 
adatoms and its neighboring carbons can be written as 
follows :" 1 

tfc = -£ ^2 t\m\CL{Rl)dm(Rl) + h.C., (12) 

I m— 0,±1 

where ClJRj) = £j =1 ^e^- 1 )/ 3 ^.^). As to 
the hopping amplitudes, to is real and t\ = t- 1 purely 
imaginary; the index / labels the plaquettes where the 
adatoms sit on the H position. Finally, the Hamiltonian 
describing the adatoms reads: 

H a = J2dHRi)h a d(Ri ), (13) 


and d^(Ri) = [dj(f?/), dj(i?/), ^^(i?/)] . The param¬ 
eters e m , A so an d A ' so are associated to crystal field ef¬ 
fects and spin-orbit coupling in the adatom p-orbitals. 

After pertubatively eliminating the adatom orbitals,"' 
the low-energy Hamiltonian in the vicinity of the two 
Dirac points near K = ( 47 r/ 3 -\/ 3 ao 7 0) (where ao = 1.42 
A is the distance between two neighboring carbon atoms 
and K = —K') reads: 

H 0 (k) = hv(r z a x k x + a y k y ) + Mr z a z s z , (15) 

where the last term T z a z s z is the so-called Kane-Mele 
(or intrinsic) spin-orbit coupling. The Fermi velocity is 
v = 3tao/(2h) and the band gap M = 3-\/3A so (being 
t = 2.7 eV in single-layer graphene). A linear Rashba 
correction to Ho(k) of the form 

R-Cr = 0R{k X Sy kyS X ) T (5r \k X [S X (JyT Z Sy(7 X ) 

-k y (s x a x + s y a y T z )] H- (16) 

with (3r = 9A^ao/2 is also induced by the adatoms. 
However, this correction is not important near the Dirac 
points at K and K' and can dropped because it yields 
subdominant corrections to the band structure, which 
otherwise realizes the Kane-Mele model. 11 For Indium 
coverage, A so ~ 8 meV and Xr = 20 meV, so M = 42 
meV and Ar/oo = 90 meV. The above continuum model 
provides an accurate description of the electronic struc¬ 
ture for energies well below the cut-off E c = Hv/ao — 4.1 
eV. 


IV. PHONON-SO COUPLING 
A. Group Theoretic Analysis 

In order to obtain the form of the phonon-SO interac¬ 
tion, the phonon displacement fields must be combined 
with electron spin and orbital operators in order to yield 
a scalar under the symmetry group Cq v . Furthermore, 
symmetry requires that the phonon field and the elec¬ 
tron operators transform according to the same irrep. In 
this work, we are interested in the coupling of the n- 
electrons near K(K '), which transform according to the 
8 -dimensional irreducible representation of C” v . We take 
three different sets of 2 x 2 Pauli matrices a, r and s to 
define the physical sublattice isospin, valley pseudospin, 
and real spin degrees of freedom respectively and these 
matrices act on eight-component spinor: 

/ c A ,K,t(k) \ 
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TABLE II: Classification of possible SO coupling terms under 
the C^ v . All these terms preserve time-reversal symmetry 


The possible SO operators transforming according to ir- 
reps of Cq v that are even under time-reversal symmetry 
are displayed in table II. 24 ’ 2 '’ 28 Note that the operators 
that contain s z are even under mirror symmetry opera¬ 
tion, z -A —z with respect to the graphene plane. On 
the other hand, the operators containing s x or s y are 
odd. The different forms of these operators for different 
choices of the spinor basis are discussed in Refs. 24,27,28, 
together with the details of the derivation of the table re¬ 
sults. 

The phonon SO coupling can be expanded in powers 
of the phonon displacement fields and its derivatives. In 
what follows, we shall focus in the long wavelength limit, 
that is, we assume that the momentum transferred by the 
phonon (q) is small (i.e. \q\ -C a^ 1 ). In what follows, 
we consider only interactions involving one phonon and 
derivative couplings will be neglected. In this regard, we 
recall that, at the T point, a uniform translation of the 
acoustic (Ei ) phono does not affect the electron dynamics 
and therefore, the coupling to these phonon modes is 
necessarily derivative and can be neglected for \q\ 0. 

In addition, we note that there is no term in Table II 
which transforms according to the E' 2 and respects the 
time reversal symmetry simultaneously. Therefore, the 
leading term from E' 2 must be dependent quadratically 
on the atomic displacement and it will be neglected in 
out treatment. 

Using the above results, we can classify the electron- 
spin to phonon coupling into two types: inter- and intra- 
valley couplings. By combining the phonon displacement 
field operators with the SO operators from table II that 
transform according to the same irrep, we can obtain 
the leading-order couplings between the phonons to the 
electron spin. Thus, the inter-valley coupling takes the 
form: 

HpSO,K = A K CFy(—UK,aSyTy — U K ,bS x Ty 

T 'U , K, c SyT x 4“ HK^dSxTx) ■ 


Whereas the intra-valley coupling reads: 

Hpso, r = Cr(ur,is x — ur,iiSy)a z T z 

+ Dr(+ur,ia y — ur,na x r z )s z . (19) 

In the above expressions, we have restricted our atten¬ 
tion to the lowest frequency optical mode transforming 
according to the G' irrep from K(K') and the E± mode 
from T, ignoring higher frequency modes (such like the 
E[ or the other G' mode). Although the coupling of 
between electron-spin and those high-frequency modes 
is not expected to vanish, their relatively high energy 
yields an scattering rate with the electrons that is (more 
strongly) exponentially suppressed at low temperatures. 
In other words, those phonon modes are less important 
in the low-energy sector that concerns us here and, there¬ 
fore, we shall not discuss them any further. 

As to the optical mode transforming according to E' 2 
near the K(K') point, even though its frequency is small 
compared to other optical modes, we shall not consider it. 
The reason is that the lowest order coupling between the 
electron-spin and this mode does not preserve time rever¬ 
sal symmetry. This can be seen from table II: It is not 
possible to write down a time-reversal invariant operator 
involving the electron spin that transforms according to 
the E 2 irrep. 

B. Microscopic derivation 

In what follows we shall extend the tight binding model 
introduced in Sect. Ill beyond the rigid-lattice (Born- 
Oppenheimer) approximation. This will allow for an al¬ 
ternative (re-)derivation of the form of the electron-spin 
to phonon couplings. In addition, we shall be able to es¬ 
timate the strength of couplings of Eqs. (18,19), which is 
not possible from the symmetry analysis provided in the 
previous section. 

In order to obtain the electron-phonon coupling, we 
need to expand the hopping amplitudes in Eq. (10) as a 
function of the atomic positions. Thus, we obtain two 
types of electron-phonon couplings. The first type re¬ 
sults from the modulation in length of the carbon-carbon 
bonds and is also present in pristine graphene. It can be 
written as follows: 18 ’ 26 



where l is the interatomic distance and b rir2 = b r5ri = 

( i 2^0: ^7.3 r 2 — — (0, l)uo nnd — 

b rir6 = i)ao are the vectors connecting a car¬ 

bon atom on the B sublattice to its nearest neighbors 
on the A. The second type stems from the modulation 
of distance between the adatoms and their neighboring 
carbons, and it is described by a term of the form: 

+ 1 Qy. 

H p c h = -J2 Y ■ Au rj CURr)dm(Ri)+ h.c., (21) 

I m= — 1 


(18) 
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where S ri = -S r4 = (0, —l)a 0 , 5 ra = -6 re = (^,^)a 0 

and 5 rs = — 5 r2 = (—-^,|)ao are the vectors con¬ 
necting the adatom to its neighboring carbons, and 
A u rj = u rj — u o with u rj and Uq associated to the dis¬ 
placement operator of carbons and adatoms respectively. 
Both types of electron-phonon couplings contribute to 
the electron-phonon scattering rate but only the latter 
one gives rise to the spin-phonon coupling in which we 
are interested here. 

In order to study the coupling within the low en¬ 
ergy sector, we can integrate out the adatom’s degrees 
of freedom by performing an unitary transformation ’ 
(see Appendix B). Hence, the strength of the phonon 
to electron-spin coupling can be obtained to be (see 
Eqs. (18) and (19)) 


and 


A k 


C r 


D r 


3(2 + a 2 )A( Q |t 1 ff| 
2vV m + "ia/2ae 0 ei 
3(2 + a 2 )A so \ti^f-\ 

4 y/m + m a /2ae\ 


( 22 ) 

(23) 


{Sab' - 3<WPT4 + 6)A( 0 |tiff | 
V2m a eoei(fo' — \Jb 12 + 4)\/ 4 + b' 2 + b’\/b' 2 + 4 


(24) 


where b' = b/a = — A + ^ k+ ^, with a = yjm a /m. 

Let us next estimate numerically the magnitude of 
these couplings. In pristine graphene, f3 = dlnt/dlnl ~ 
2 to 3. 18 Thus, we choose /3 = 2 and estimate dto/dl ~ 
(3t 0 /ao and dt\/dl ~ /3ti/ao- Using the parame¬ 
ters given in Ref. [16], we estimate the intra-valley s z - 
conserving electron-phonon coupling to be Dr ~ 46 
meV/A \fm , intra-valley spin-flip electron-phonon cou¬ 
pling Cr ~ 98 meV /A \fm, and the inter-valley spin- 
flip coupling Ak ~ 28 meV/A ^fm. In addition, using 
Eq. (18) and Eq. (19), we can study the effective electron- 
electron interaction by integrating out the phonon modes 
ur,i,UT.n,UaiUb,u c ,Ud- The dominative effective inter¬ 
actions induced by phonon are: G 2 /2w 2 ~ 2.6 meV 
_Dp/2u; 2 ~ 0.57 meV and A 2 k /2uj 2 k ~ 1.3 meV. 


V. APPLICATION: QUASI-PARTICLE 
SPIN-FLIP RATE DUE TO PHONONS 

A. Formalism 

The electron-spin to phonon coupling displayed in Eqs. 
(18) and (19) is the main result of this work. It can be 
applied to different problems, such like the calculation of 
spin transport in the system as well as the study of su¬ 
perconducting pairing instabilities in ad-atom decorated 
graphene. However, in this section, we focus on comput¬ 
ing the spin-flip rate of electrons and holes. This will 
provide a quantitative estimation of the importance of 
this coupling as a function of the quasi-particle energy 
and system temperature. 


Indeed, the existence of the spin-electron phonon cou¬ 
pling implies that phonon emission and absorption can 
flip the spin of electronic quasiparticles. The spin-flip 
rate, T, can be obtained from the imaginary part of 
the electron self-energy (E) as T = —Im E. The self¬ 
energy is related to the single-particle Green’s function 
(G) by means of Dyson’s equation: G _1 = G^ 1 — E, 
where Go is the non-interacting electron Green’s func¬ 
tion. The quasiparticle decay rate, T, is a sum of dif¬ 
ferent contributions T = T ee + Tdi s + T e _ p h + • • •, which 
contains contributions from electron-electron interactions 
T ee , disorder Tdi s , lattice vibrations T e _ p /j, and other in¬ 
elastic mechanisms. Here we focus on r e _ p h term, which 
contains contributions from two different channels: the 
spin-flip channel (T fllp ) and non-spin-flip channel (T non ). 
The non-spin-flip channel contribution arises from the 
electron-phonon couplings of various kinds. On the other 
hand, the spin-flip contribution is exclusively generated 
by the spin-phonon coupling that has been discussed 
above. Note that such phonon-mediated spin-flip pro¬ 
cesses are different from those due to electron-electron 
interactions and (magnetic/SOC) impurity scattering be¬ 
cause they typically have a much more pronounced tem¬ 
perature and energy dependence due to the optical nature 
of the phonons involved. Therefore, the understanding of 
this spin-flip processes becomes necessary when thinking 
about possible device applications, which must be oper¬ 
ated in a finite temperature. 

In order to investigate the spin-flip rate due to 
phonons, we first compute the (non-interacting) quasi¬ 
particle dispersion by diagonalizing the Hamiltonian in 
Eq. (15). This can be achieved by means of a unitary 
transformation, U (see Appendix C for the details): 


tf(fc) = {7 -1 $(fe) 


/ c +tK ',i(k) \ 

C+,K,l(ty 

c+,x',t(fc) 
c+,k, t(fc) 
C- t K',i{k) ’ 

C-,K,±{k) 

C-,K>, t( fc ) 

V c_ K,t(*0 / 


(25) 


where $(fc) is the original basis given by Eq. (17). As a 
result, the unperturbed Hamiltonian becomes diagonal, 

i.e. 


U ] H 0 (k)U - nN 


£k, + 1 4x4 o 

0 £fc,_/4x4 


N 


(26) 


where £k.± = ±\/M 2 + H 2 v 2 k 2 — fx. Here fi is the chem¬ 
ical potential and l 4 X 4 is the 4x4 unit matrix. 

Since the electron spin is coupled to the phonon by 
means of (18) and (19), which are both linear in atomic 
displacement (i.e. phonon field), the decay rate including 
spin-flip and non spin-flip processes can be obtained by 
calculating the imaginary part of the electron self-energy 
to lowest order in the electron-phonon interaction . We 
shall neglect the renormalization of the quasi-particle 



















energy, which effectively renders our approach entirely 
equivalent to a Fermi-Golden rule calculation. Thus, the 
decay rate for a quasi-particle (or quasi-hole) with energy 
lo (referred to Fermi energy) and band-index ± (where 
+ stands for conduction and — for valence band), valley 
K(K') and spin f at momentum k is given by 

r h ,» = 

p,r]i a 

j^(l “b n(w Q ) f (^fc—p,?^)) & (tU (^fc— p.r/i T" ^a)) 

-t-(u(ce a ) T f (^fc—p.r^i)) 3(yi (^fc— p,T]i ^a)) ■ 

(27) 

where ry = (±, A"(AT'), t (4-)) denotes the set of quasi¬ 
particle quantum numbers (band, valley, spin), and N is 
the total number of unit cells. L a stands for the 4x4 
matrices in Eq. (18) and the two terms in Eq. (19) whose 
couplings are J a = Ax, Cp, £>r- The phonon frequencies 
uj a = uJr,ux are given in Eqs. (4,7); /(w) and n(w) are 
the Fermi-Dirac and Bose-Einstein distributions, respec¬ 
tively. 



FIG. 4: Quasi-particle spin-flip rate due to the electron spin- 
phonon coupling, as a function of quasi-particle energy cj 
when T = 10A, the band gap is A = 0.1 eV and chemi¬ 
cal potential around 0.11 eV, measured from the middle of 
the spin gap. The blue line is for the total spin-flip rate 
r fllp (ta) = rf4 P ra (o;)-|-rf^ p er (a;), red dashed line for intra-valley 
spin-flip rfx a ( to) and green dash-dotted line for intervalley 
spin-flip rfX r H. The inset shows the band structure near 
a Dirac point with a spin gap with a black thick line asso¬ 
ciated with the location of the chemical potential. For this 
plot, we use the parameters of Indium coverage discussed in 
context. Note that r® lp (w),rf at p raiA (a;),rf at p er!A (w) can have 
values only when \(oj + /i) > M, but for convenience we plot 
them together to show overall features among the whole range 
of u>. 


(a) 




tu(eV) 


FIG. 5: Spin-flip rate F filp (u;) as a function of the quasi¬ 
particle energy uj. (a) for chemical potentials /r = 0.11 eV 
and 0 eV for T = 10A". (b) for various temperature T = 
10, 50,100A for yr = 0.11 eV. The values of the remaining 
parameters are the same as in Fig. 4 


B. Quasi-particle Spin-flip rate 

The spin-flip rate of a spin polarized quasi-particle, 
say, an electron injected at any given valley with spin 
up and energy hu (relative to the Fermi energy), may 
be easier to determine experimentally. Note that, as 
mentioned above, we are interested in the contribution 
from the spin-flip rate as obtained from the electron spin- 
phonon coupling derived above (processes that do not 
flip the spin can also be induced by other mechanisms 
such like the normal electron-phonon coupling). Thus, 
we shall retain the contribution from the following four 
matrix elements: (WL a U)( a -)-) (b 4 .) in Eq. (27), where 
a,b = (±, K/K'). 

For the current example, the band structures are 
isotropic for momentum k and symmetric between un¬ 
like spins at low energy. Therefore, after some manip¬ 
ulations (see Appendix D for details), the spin-flip rate 
can be shown to be independent on the direction of mo¬ 
mentum k and the spin of the quasi-particle. As a con¬ 
sequence, the spin-flip rate of a spin polarized quasi- 
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particle around Fermi surface can be written as 

rj'» = r?i» + 

27r 

= — [2A^L\(uj,lok ) + C 2 L\(u>, wr)] ,(^ 8 ) 

where A = ± refers to the band index. r?^ er A and r?^ ra A 
are the contributions of the inter-valley and intra-valley 
phonon SO couplings (see Eqs. (18) and (19)) respec¬ 
tively. Here L\(u>,uj a ) is defined by (see Appendix D) 


L\(u>,u>a) = Q(Xhco' — M ) 


n 


2Nhu} a Trv 2 


r 

1 

i . , M 2 

x { 

1 +n(w Q ) - /(w - w a ) 

n\u u a \ 

1 1 h\u)\i 


x e(*v -u a y - m 2 ) 

-f* TliJVjJa) T f (cu T u^a) h^to T uj a | — 
xQ(H 2 (uj' + w Q ) 2 — M 2 )} , 


M 2 

h\ui'\. 


(29) 


where H is the system size (area) and N is the number 
of unit cells. ©(*) is Heaviside step function, and u/ = 
ui + )jl. The ratio between coefficients in front of A 2 K and 
Op in Eq. (28) result from the fact that Eq. (18) contains 
four terms proportional to s x ,s y , which can cause spin- 
flip, whereas Eq. (19) contains only two such terms. 

In Fig. 4, we have plotted the quasi-particle spin-flip 
rate as a function of the quasi-particle energy, ui (rela¬ 
tive to the Fermi energy). The Fermi energy (chemical 
potential at T = 0), /z, is measured from the middle of 
the band gap. We have set the temperature T = 10 K. 
The band structure near a Dirac point is shown in the in¬ 
set. The quasiparticle spin-flip rate vanishes between ui — 
—0.16 eV and w — —0.07 eV because quasiparticle propa¬ 
gation is not possible within the band gap. Furthermore, 
energy conservation suppresses the phonon emission at 
low temperatures: In order to emit a phonon, the elec¬ 
tron quasi-particle energy (ui > 0) must satisfy the con¬ 
dition that w = £fc = £fc_q +w^(r). For a quasi-hole (i.e. 
for ui < 0) we must have ui = —£k = — £k-q + ^k/t- 
In addition, Fermi statistics further prevents an electron 
from emitting a phonon if the final state is below the 
Fermi energy. As a consequence, since the electron spin- 
flip is caused by optical phonons with energies u>k and 
wr > ujk (for q —> 0), the quasi-particle spin-flip rate ex¬ 
hibits a two-step variation close to w = 0 in Fig. 4. The 
phenomena applies to holes: Because of Fermi statistics, 
a double step variation occurs only near the top of the 
valance band and near (but below) the Fermi energy. The 
variation of the spin-flip rate with the chemical potential, 
/i, is shown in In Fig. 5(a) for T = 10 K. The temper¬ 
ature dependence, however, is displayed in Fig. 5(b) for 
(i ~ O.lleH. The double step structure still occurs for 
different values of /j,. However, increasing the tempera¬ 
ture smears out this structure. 

Finally, we briefly remark the effects of the Rashba 
coupling (cf. Eq. (16)), which has been neglected when 


calculating the spin flip rate due to the coupling with 
phonons. There are two different influences of this term 
on the spin motion: First, the presence of the Rashba 
term in the band structure leads corrections to the in¬ 
elastic spin flip rate through the change of the density 
of states (i.e. the slope of the energy dispersion at 
low energies). The effect can be estimated by replac¬ 
ing the prefactor 1/v 2 in Eq. (29) by 1 /(v ± Pr/H ) 2 « 
1/v 2 2/3^/fi 2 ri 2 . Thus we see that the correction arising 

from the linear Rashba term is ~ fix/hi) 2 w 10 3 , which 
can be safely neglected. 

The second effect results from the fact that the elec¬ 
tron spin is no longer a good quantum number when the 
Rashba term is included. In other words, an injected 
electron whose spin is pointing up in the out of the plane 
direction, for instance, will naturally precess even without 
coupling to the phonons. In other words, the spin state 
will experience Rabi-like oscillations. Nevertheless, this 
Rashba-induced precession does not lead to energy loss 
of the injected electron and therefore, it should not af¬ 
fect the inelastic spin flip rate mediated by phonons that 
we have computed above. We can estimate the preces¬ 
sion rate (T/j) by diagonalizing the electron Hamiltonian 
including the order correction to order /3r/Hv: 


Fr/M 


1 /3r A* 
hhv M 1 


(30) 


where fx is the chemical potential. If we use the same 
parameters as in Fig. 4 taking /r ~ 100 meV, we find 
r# ~ 10 4 ns -1 , which is much larger than the spin 
flip rate (Tfiip) due to spin-phonon interaction. There¬ 
fore, in an experiment, the spin orientation of the in¬ 
jected electron will precess very fast due to the spin- 
orbital (Rashba) coupling. However, such precession will 
be damped by inelastic processes arising from e.g. the 
coupling to the phonons. Thus, the double step struc¬ 
ture as a function of the quasi-particle energy that we 
found above should still appear even in the presence of a 
Rashba-type spin-orbit coupling. 


VI. DISCUSSION 

We have obtained the electron-spin to phonon cou¬ 
pling for a single-layer of graphene uniformly covered 
with heavy adatoms. The coupling is obtained using 
both symmetry arguments and a simple tight-binding 
model. However, when considering the experimental con¬ 
sequences of our calculations, a major concern is that, as 
it has been found (see e.g. Ref. [19,31]), adatoms tend 
to cluster rather than uniformly covering the graphene 
layer. Indeed, it has been discussed by a number of au¬ 
thors recently, some of the features of the (topological) 
band gap may be still observable in the quasi-particle 
spectrum or in the transport properties ’ provided 
clustering is not very strong . In such scenario, we ex¬ 
pect that, even in the presence of a sizable randomness 
in adatoms coverage, a coupling between the electron 
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spin and the lattice dynamics exists. Assuming that a 
large faction of the adatoms are located in the H po¬ 
sition, we can expect that the randomness will lead to 
some sizable broadening of the phonon spectrum (espe¬ 
cially of the shorter wave-length phonons at the K and 
K' points). The latter will result in a rounding up of 
the sharp features of the energy dependence of spin-flip 
rate shown in Fig. 4 and 5, even at zero temperature. 
For a clean and uniformly covered single-layer graphene 
system, the existence of these sharp features and the ex¬ 
ponential suppression of the spin-flip processes that they 
imply for quasi-particles near the Fermi level seem to im¬ 
ply that the edge transport may be essentially ballistic 
temperatures, if disorder exists in the adatom coverage, 
this may not be longer the case. The quantification of 
this effect is currently under investigation . 

In addition, our results also open the door to research 
into superconductivity in decorated graphene. Indeed, in 
recent time, studies of layer materials such as graphene 
heterostructures separated by hexagonal boron nitride 
and molybdenum disulfide were proposed, when the 
carrier concentration is high and the screening of long 
range Coulomb potential is strong, to have a supercon¬ 
ducting phase, at which the cooper pairing has opposite 
signs in different valleys. It becomes more interesting in 
the graphene decorated with heavy adatoms because the 
couplings between the electron-spin and lattice vibration 
can provide some more unusual effective spin-spin inter¬ 
actions. Competition between electronic repulsion and 
the effective spin-spin interaction with tunable charge 
density can offer a new platform to show rich unconven¬ 
tional superconductivity phases. 


VII. SUMMARY 


In summary, we have used symmetry arguments and 
tight-binding models to obtain the coupling between the 
lattice dynamics and the electron spin in adatom dec¬ 
orated graphene. Since the vibration mode transform¬ 
ing according to the G'-representation at K and K' and 
the mode transforming according to E\ at T are the 
lowest energy phonons coupling to the electron-spin, we 
have focused on them. The strength of the phonon to 
electron-spin coupling constants have been estimated us¬ 
ing a tight-binding model. Finally, in order to study 
the effect of this coupling on the dynamics of electronic 
quasi-particle, we have computed the phonon spin-flip 
rate. Our results could be relevant for the understand¬ 
ing of the temperature dependence of spin-transport in 
graphene decorated with heavy adatoms. In particular, 
the spin-flip rate may be measurable by performing spin- 
polarize scanning tunneling spectroscopy of the system. 
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Appendix A: Dynamics matrix of the center force 
model 


In this appendix, we provide the details of the deriva¬ 
tion of the central-force model introduced in section IIB 
(cf. Eq. (1)). By working in Fourier space, we obtain the 
phonon eigenfrequencies and eigenmodes of the model. 
Upon inserting Ui(R) = J2 P e lp ' ( ' R+<5 ^ } Ui(p) into 
Eq. (1), the elastic energy term, J2 t j •( M *-'iu)) 2 , 

becomes: 




-A 

0 

G* 

F* 

C* 

G*- 

0 

A 

F* 

D\ 

G* 

d* 2 

c\ 

Ei 

A 

0 

Ef 

G 

Fi 

D 1 

0 

A 

G 

g 2 

C 2 

G 

f 2 

G* 

B 

0 

G 

d 2 

G* 

g 2 

0 

B 


(Al) 


where 


u T (q) = (. x A (q),yA{q),x B (q),yB(q),x 0 (q),yo(q )), 




(A2) 

A = 

3 3,, 

2 k + 2 k ’ 

(A3) 

B = 

3 k', 

(A4) 

C\ = 

-^k[e ipbl2 + e ipb a 4 ], 

(A5) 

c 2 = 

~kt[e ipmSa +e ip -**\, 

(A6) 

Di = 

_I k[e ipbl2 + e lpb3i + 4e lpb32 ], 

(A7) 

d 2 = 

--k'\e ip - SB + e ipSa + Ae ipSl ], 

4 

(AS) 

F\ = 

~^k[e ip ' bl2 - e ip ' b34 ], 

(A9) 

f 2 = 

--F[e tp - &2 + e lp - s,i ] = G 2 *, 

(A10) 

G = 

k'[e ipS3 — e ip ' S5 }, 

(All) 

G 2 = 

D*. 

(A12) 


Here <5i = -<5 4 = a o (0,-l),5 3 = -S 6 = ao(-^, |) 

and < 5,5 = —S 2 = ao(—^, ^) are the vectors joining the 
adatom to its neighboring carbon atoms; b 12 = b 54 = 
a o(-^!- 5 )) b 32 = b 56 = a 0 ( 0 , 1 ) and 634 = b 16 = 




ao(-^, — ^) are the vectors joining a carbon atom on the 
B sublattice to the neighboring carbons on the A sublat¬ 
tice (see Fig. 1). 


Appendix B: Perturbation theory 


In this Appendix, we outline the details of the per¬ 
turbation theory ' 1 ’ 1 ' 1 that allows us to integrate out 
the adatom orbitals and obtain the proximity-induced 
electron-spin to phonon coupling. The starting electronic 
Hamiltonian matrix consists of the minimal tight-binding 
model, Eq. (10), discussed in section III and the electron- 
phonon coupling, Eq. (20) and Eq. (21), in section IVB. 
By means of a Fourier transformation, 

dm(Ri) = -7^J2 e lk Rl dm(k), (Bl) 

CrARl ) = ^E eifc ' (Rl+S)c ( fe )’ ( B2 ) 

u rj (Rj) = ^E eifc ' (R/+S M fc )> (B3) 

v k 


the electronic Hamiltonian matrix in the momentum 
space can be written as follows: 

H = ^2i>iH(k,k-p)ip k -p, (B4) 

k,p 


where 


H(k,k — p) 


hg(k, k — p) 

T(k,k-p) 


T{k,k-p) 
h„(k, k — p) 


, (B5) 


and 


V'fc = [c A {k),c B {k),di{k),d 0 (k),d- 1 (k)} . (B 6 ) 


Here h g (k,k — p) is the Hamiltonian matrix just coming 
from graphene (see Eqs.(ll) and (20)), but its explicit 
form is not important for the calculation sketched below. 
h a (k, k — p) represents the Hamiltonian of adatoms: 


h a (k, k-p) 


ei -f- A sos z \/2A' s 0 s 0 

V^A' so s + eo %/2A 'so s 
0 %/2As 0 s + £ i — A sos" _ 

(B7) 


T(k,k — p) and T(k,k — p) come from, Eq. (12) and 
Eq. (21) respectively, which allow electrons to tunnel be¬ 
tween the adatoms and its neighboring carbons with the 
influence of adatom’s displacement. After some straight¬ 
forward algebra, we have T(k, k — p) = To(k)5k,k-p + 
T u (k, k—p), and T(k, k—p) = T^(k)5 k ,k- P +T u {k, k-p), 
where T u (k, k — p) can be obtained by 


T u (k, k — p) = Tt(-p, k - p). 


11 


Here we use p —► —p first and then k —> k p on 
T^(fe, k — p). T 0 (k) is defined as: 


T 0 (k) = 


Vi(fc) V 0 (le) y_i(fe)' 
U\{k) U 0 (k) J7_i(fe) 


(B9) 

(BIO) 


with 


Vm(k) — t\m\ ^ ^ 


0 i{n— l)7r/3 iS n k 


n= 1,3,5 

U m (k) = H! m (fc), 
and at the same time 


(Bll) 

(B12) 


T u (k, k- p) = 

Vi A(k,k-p) V 0 A(k, k-p) V-iA(k, k-p) 
Vi B (k,k-p) V 0B (k. k-p) V-i B (k,k-p) 


(B13) 


where 

VmA{k, k- p) = 


1 ^l m l [ xa{p ) [ v ^ c i ( 27| 7 3 ) m c ~ ilS 3 ( fc -p) 
y/6 N dl 2 


_ }/3. e i(4,n/3)m £ -iS 5 -(k~p)^ | x„(p)[ c i ( 2n / 3 ') m c ~ iSs ' k 

+ ^gi(4vr/3 + ^ ( p ) ’( fe “P) 


i(2ir/3)m ^-iS 3 (k-p) _j_ £ e i(4ir/3)m g-i *5 (h-p) j 


+ 2 6 

/ \ r i —iS-\-k 1 i(2ir/3)m -iS 3 -k 1 i(47r/3)m — idc-Zen 

+ yo(p)[+e --e y ’ ’ e 3 --e y 1 ’ e ]} 


and 


VmB{k, k- p) = 


(B14) 


1 dt\m\ j( 7r /3) rn ~-iS 2 .(k-p) 

Vm dl iXB[p)l 2 


c H5ir/3)m c ~ig 6 .(fe-p)j ^ (p) [_ e i (.^/3)m e ~i&2 ■ *» 

'/S i(5Tr/3)m -iS 6 k-i , / \ [_ 1 +i(n/3)m -iS 2 -(k-p) 


+ — 6 


"e ■ ] + ys(p)[--e 


_j_ giTrm e -iS 4 -(k-p) _ /3)m ^-iSQ-{k-p )j 

/ \ r i 1 i(7r/3)m -iSo-k iirm —iSa-k 

+ Vo(p)[+ 2 e e 2 -e e 

_j_ lgi( 5 ’ r /3)m e -i«5 6 -fej j 


(B15) 


To integrate out adatom’s degrees of freedom, we can 
choose a proper transformation H —► H = e~ s He s to 
make h g and h a decoupled. When 


S = 


0 M 
-Aft 0 


(B16) 


the Hamiltonian matrix is transformed to 


H = e~ s He s 


H 


n 


Ml/i, + T - h a M^ 


hgM + T - Mh a 

H 22 


- k,k-p 

(B17) 


(B 8 ) 
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The off-diagonal terms can be eliminated when requiring 

f = —M^hg + h a M\ (B18) 

T = Mh a - h g M , (B19) 

or equivalently 

Aft = h~ lr T + h~ 1 M+ g , (B20) 

M = Th~ 1 +h g Mh~ 1 . (B21) 

This transformation yields corrections to h gi such that 

h g = h g - TM t - MT - i h g MM t - ^MM^h g 

+ Mh a M f + 0(M 3 ) (B22) 

~ hg-Th-'T, (B23) 

where we employ i(— h g MM t + Mh a M t) = i TM t and 

+ Mh a M f) = ±MT. Following the as¬ 


sumption used in Ref. [16] that the spin-orbit splitting 
of adatom orbitals is small compared to the crystal field 
effects, we can keep only the leading order terms in 
A go and A so and make h~ l ~ h^ 1 — h^ 1 hih^ 1 , where 
h a = h o + hi and hg = diag(ei, eo, ei). Thus, to the lead¬ 
ing order of phonon displacement fields we obtain 

h g ( k , k S pi n ~ To ( hi h q T u (/c , k p} 

+ T u (k,k - p)h~ 1 hih~ 1 T^(k - p). 

(B24) 

Here only terms relevant to electron’s spin are shown. 
Eq. (B24) can help us to compute strength of the phonon- 
SO coupling of Eqs. (18,19). 

Appendix C: Diagonalization of the non-interacting 
electron Hamiltonian 


I 

In this appendix, we will show the electronic Hamiltonian given in Eq. (15), 

H 0 (k) = v F h(r z a x k x + u y k y ) + Mr z a z s z , (Cl) 


in the basis, 


$(fc) 


/ c a ,k. t(fc)\ 
cb.k, t(fc) 

CA, K'+k) 

CB, K’, t(fe) 
CA,K,i(k) ’ 

cb,k, i(k) 

CA.K', J.(fc) 
\CB,K',i{k)J 


(C2) 


can be diagonalized by an unitary transformation E{k) = U 1 Ho(k)U, where the unitary matrix is 


' 0 

0 

0 

M+R 

F + 

0 

0 

0 

M-R -i 
F_ 

0 

0 

0 

hvke i& k 

F + 

0 

0 

0 

hvke id k 

F_ 

0 

0 

M-R 

F_ 

0 

0 

0 

M+R 

F + 

0 

0 

0 

hvke~ iG k 

F_ 

0 

0 

0 

hvke 

*+ 

0 

0 

-M+R 

F_ 

0 

0 

0 

-M-R 

F + a 

0 

0 

0 

hvke^k 

F_ 

0 

0 

0 

hvke^ks 

F + 

0 

0 

M+R 

F + 

0 

0 

0 

M-R 

F_ 

0 

0 

0 

hvke 

L F + 

0 

0 

0 

hvke~ l ®k. 

F_ 

0 

0 

0 


(C3) 


with tan0fc = k y /k x , R(k ) = yj A T 2 + h 2 v 2 (k x + k y ), 
F + (k) = yj (M + R(k )) 2 + h 2 v 2 {k x + k y ) and 

F_(k) = — R(k)) 2 + h 2 v 2 {kl + k%). After this 

unitary transformation, the electronic Hamiltonian in 
the grand canonical ensemble becomes a diagonal energy 


eigenvalue matrix, 


E(k)~ p 


( R(k ) - p)hx4 0 

0 (-R(k)~n)h x 4 

£fc, + l4x4 0 

0 Ek, — 1 4x4 ’ 


(C4) 

(C5) 


where £k,± = M 2 + h?v 2 k 2 — p. Here /r is the chem- 
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terms: (J 7 t £“[/)( Q ,t) > (& 4 ) from Eq. (27), where a,b = 
(±,K/K'). Therefore the spin-flip rate of a quasiparti¬ 
cle with rj = ().)) quantum numbers can be 

written as 

To compute the spin-flip rate, we can straight¬ 
forwardly select contribution from the following four 

I 


ical potential and / 4 X 4 is the 4x4 unit matrix. 

Appendix D: The spin-flip rate 


pflip _ 2-7T y- 2 A 2 k 
lfe .± - h 2^ 


h ' 2 ujk 
v 


g+{u K ,k,p ) +g-(uj K ,k,p) ± 


-g+(u K ,k,p) + g-(u K ,k,p) 
R(k)R(k — p) 


M 2 - -k\k - p\v 2 cos( 6 » fe - d h -p) 


+ jfl+^r, fc, p) + g-(ur, k, p) ± ff+(^fc^ + g ^ r ’ fc,P ^ [.M 2 + fc|fc - p\v 2 cos (8 k - 6 k -p)\ j , (Dl) 


where 


1 1 

g±(oj a ,k,p) = — [l + n(« a ) -/(efe- Pl ±)]£(w - (efe_p,±+« a )) + — [n(w a )+ /(efe_p,±)]<Hw-(efe- P ,± - Wet)), (D 2 ) 


N 


and R(k) = yjM 2 + H 2 v 2 {k 2 + k 2 ). Because the band 
structures are symmetric between unlike spins and unlike 
valleys at low energy, the spin-flip rate from the quasi¬ 
particle of channel y is dependent on the energy of quasi¬ 
particles, ui = £k.±, but independent on the spin and 


valley labels. In the following calculation we would like 
to deal with the summation on the internal momentum p. 
By rewriting the summation as integral over the internal 
momentum p, Eq. (Dl) becomes 


pflip _ 2 tt 

ifc ’± ~ h 


dP l rfe 2 w^‘ P { ffA<WK ’ fc ’ P) 


Xgx(uJK,k,p) 


M 2 — —kv 2 (k — p cos 9) 


R(k)R(k — p) 

+ 2 iHj£yf 0 d pf o d e^p{g^r,k,p)T ll^L P) p) [M 2 + kv 2 (k-pco S 9)\ 


(D3) 


2 tt 

T 


2A 2 k V r dp [ \y- 


1 -p 


2 Uk yj\ — y 

2tt^ 2 v- 2 d r°° , r 1 , 1 -p 
+ T Cr S w/o *1 dv ^ 


2 wr x/l — y 2 


gx(oJ K ,k,p )-F 
gx{u T ,k,p) =f 

I- 


\gx(uj K ,k,p) 


R(k)R(k — p) 
Xgx(uJr,k,p) r 2 , ,...2 


M 2 — -kv 2 (k — py) 


R(k)R(k — p) 


[M 2 + kv 2 {k — py)] 


(D4) 


where cos( 0 *, — 0 fc_ p ) = (k—pcosd)/\k—p\, y = cos 9 and where u/ = ui + p and 

—h 2 (u}' =p uj a ) 2 + M 2 + h 2 v 2 (k 2 + p 2 ) 


the azimuthal angle part is replaced as d9 = — 7 ==dy. 

v 1— y 2 


Let us introduce the function: 

e k - P ,\(y) = Ai?(fe — p) 


J/=F,a 


2 h 2 v 2 kp 


(D7) 


= XxJ M 2 + h 2 v 2 {k 2 +p 2 — 2 kpy), (D5) 

with A = ±. Thus, the delta function of Eq. (D2) be¬ 
comes 

c/ / f€-k—p,\(y) I \\ |efc — p, A (y=F,a) I c/ \ 

5(w - ( ^r~ ± Wq)) = hv*kp 5{y - y ^ 

X 0(A h{u' =F Wa) - M) (D6) 


Furthermore the interval of the integral of p is limited 
by —1 < y^z, a < 1 (remember y = cos0). From the 
inequality, j/ TiCe < 1 , we obtain 


Pi,T,« < A < P2,^,o 


where 


Pl,=F,a = k- 


h 2 (u}' yu a ) 2 -M 2 


v 2 h 2 


(D 8 ) 


(D9) 
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and 


and 


P2,^,a — k + 


S 2 (w'Tw„) 2 - M 2 


y 2 fi , 2 


(DIO) 


with ^/h 2 {uj' =F w a ) 2 — Af 2 being real. From the other 
inequality, — 1 < y=p, a , we obtain 


P 2 , T ,a = - k - 


h 2 (u>'^Lu a ) 2 - M 2 


v 2 h 2 


(D13) 


P < P2,^Pl,^,a <P 


(Dll) 


where 


Pi,=F,« — + 


fi 2 (w' q= u„) 2 — M 2 


Therefore these inequalities limit the range of p within 


u 2 h 2 


(D12) Max.{pi lTia ,P* tTta } <P< P‘2, a- 


Bringing results of above discussions into Eq. (D4), we have (skip coefficients ^-2A 2 K and ^2 Cf- here) 


Sispr*/ dv 7t 


x=± 

2ft 


— y 2 2cu ( 


poo p — 1 

E / «»/ 

,_, ./o di 


AT(27r) 


dy 


-gx(u} a ,k,p) 
1 


(D14) 


^/1 — y2 2 hv 2 kpUJ a 


ek-p,\(y-,a) (l + n(w a ) - f(e k -p,x(y) - p))S(y - y_, a )0(Afi(u/ - uj a ) - M) 


+ \e h -p,\{y+,a) (ra(w a ) + f(ek-p,x(y) ~ p))S(y - y+, a )®(\h{u' + u a ) - M ) 


S L *{ I 1+ -“- 1 -">]-“»)-«) 


(D15) 


\u' + W a |0(l - \y+,<x\). 


and 


+ \n(huj a ) + f(XH\u}' + 0 J a \ - p)] lUJ ^ - = -0(Afi(u/ + U a ) — M) 

L J 2 tv a y/l- (y+, a ) 2 

-p 1 *9\{uK,k,p) 


(D16) 




A=± 

2D 


N(2tv) 2 


roo r — 1 

^ / dp dy- 
x=± J o -D 


— y 2 2wjc R(k)R(k — p) 
-p A 


M 2 — -kv 2 (k — py) 


(D17) 


,/rr^ 2 hv 2 k P R{k)^ { I 1 + n( ^ } - /(efe -- A(2/) - M) ] 5{y ~ ~ "*) - M ) 

+ |n(wje) + f(ek-p,\(y) - p]js(y - y + , a )0(\h(u>' +co K ) - M) 


M 2 — i fcr) 2 (fe — py) 


(D18) 


2D 


Nhv 2 kR(k){2n) 2 


A J dp j ^1 + n(w K ) - /(AR|w' - wjr| - p)J 


0(1 - \y-, a \)(2M 2 - k 2 v 2 +pkv 2 y -, a ) 


x0(Afi(w' — uk) — Af) + |^n(/kjif) + f(Xh\uj' + ujk\ — /x)j 


AlOKx/1 - ( y-,a) 2 
0(1 - |j/ + ,o:|)(2M 2 - fc 2 « 2 +pfcn 2 t/ + , n ) 

4w*V! - (y +,«) 2 


0(Afr(cj' + u>k) — M) j , 
(D19) 


r 


where \ek- P ,\(y±, a )\ = W ± w a |- In our case the opti- the integral: 
cal phonon mode uj a does not have momentum depen¬ 
dence (to the leading order), so we just need to compute 


dp 


0(1 - \y±, a \) 


L x/ 1 ~ ( y ±,<*) 2 

fP 2,±,a 

/ dp 

■/Max{p 1±ct ,p* , } 


’ p i,±,oi a/i (y±,a) 2 


rP2,±,a 

J Max{pi 
rP2,±,a 

L ?,. 


±.«>pi,±,, 


dp ^ 

} \/P 2 “ (P±,c +/3p 2 ) 2 


(D20) 

(D21) 

(D22) 


2-^* - (p±,o + fix) 2 


= 2tt/c, (D23) 
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and 


/r 2 i * 2 „, 2 u 2 


dp 


PV±,a Q(1 - \y±,a\) 

V 1 - (y±,«) 2 


rP 2 ,±„ 


'Pf,±„ 


dx 


P±,a + /?£ 


2y/x- (/ 9 ±,a + /?£)' 


(D24) 


= 27rfc 2 ,(D25) 


w h cr e w , = p±-°+^ 0 , = -ft 2 (o/±<A,) 2 +M 2 +ft 2 

wnere y±, a — p , p±, a — 2 h 2 v 2 k 

/3 = and x = p 2 . Because of Eq. (D23) and (D25), 
cos(0fc — 0/ c _ p ) in Eq. (Dl) gives no contribution to the 
spin flip rate. Therefore terms associated with uk i n 
can be simplified as 


dv 


\J\ — y 2 2 ujk 


gx(uJK,k,p) =F 


A gx{uK,k,p) 
R(k)R(k — p) 


M 2 — —kv 2 (k — py) 


ft 


2NhwK^v 2 


| [l + n(uj K ) - f{u - luk) j \h\u' - u K \ =F sgn(w' - u K ) - w^) 2 - A4 2 ) 


+ |n(fiui/r) + /(w + + ojk | =F sgn(w' + u K ) ] 0 (ft 2 (a/ + w) 2 - A/ 2 ) j 


= L±(w,wk)', 

Terms associated with air can be simplified as 


(D26) 


(D27) 

(D28) 


e ^ r d -r d «^fhM s ^ r ' kv) t mw-i) k + kAk - w) ]} ° l±(w) - tD29) 

r 


Since the energy of incoming quasi-particle is ui, we can with L±(uj,uj a ) and sgn(a/ ± lok, r) can be simplified as 
set R(k) = |o/| = \w + p\ in Eq. (D28) and (D29) and sgn(w') (M > fuxK.v in our case). Thus the spin-flip rate 
drop the k dependence of Tj.'^o;). To ensure incoming °f a spin polarized quasi-particle at a valley in a band 
quasi-particles within bands, 0(±u/ — M) is introduced denoted by A = ± is obtained to be 


r A P M = r ?nt r » + rSL, A M = ^ [2A 2 k L x (u, u k ) + C 2 L\(u, Wr )], 


-\flip 


2n i 


(D30) 


where 


L\(co,ijj a ) = 0(Afio/ — M) 


n 


2Nhuj a nv 2 


r -| 

1 . , M 2 1 

1 + n(u a ) - f(u - W«) 

n\uj u a \ 

L 1 1 n\uj'\i 


<~>{h 2 (u}' - uj a y - m 2 ) 


+ 


n(ui a ) + f(ui + w a ) 


h\u> T uo a I — 


M 2 


h\io'\ 


9{h 2 (oj' + u a ) 2 - M 2 ) 


(D31) 


In the expression of T A ip (u;), the scattering event associ¬ 
ated with u>k phonons of G' modes is identified as inter¬ 
valley spin-flip rfj er A (w), whereas the event associated 


with wr phonons of E\ modes is identified as intravalley 
spin-flip rf n ip A (w). 
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